Generation of two-photon states with arbitrary degree of entanglement via nonlinear 

crystal superlattices 
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We demonstrate a general method of engineering the joint quantum state of photon pairs produced 
in spontaneous parametric downconversion (PDC). The method makes use of a superlattice struc- 
ture of nonlinear and linear materials, in conjunction with a broadband pump, to manipulate the 
group delays of the signal and idler photons relative to the pump pulse, and realizes photon pairs de- 
scribed by a joint spectral amplitude with arbitrary degree of entanglement. This method of group 
delay engineering has the potential of synthesizing a broad range of states including factorizable 
states crucial for quantum networking and states optimized for Hong-Ou-Mandel interferometry. 
Experimental results for the latter case are presented, illustrating the principles of this approach. 

PACS numbers: 42.50.Ar, 03.67.Lx 



Quantum interference between single-photon 
wavepackets lies at the heart of novel quantum-enhanced 
communication and computation technologies. Such 
interference depends crucially on the spatio-temporal 
mode structure of the photons used. Therefore it 
is important to develop techniques that enable the 
generation of quantum states with the specific modal 
character required for different applications. In par- 
ticular, entangled photon pairs provide an important 
resource for several quantum information processing 
(QIP) protocols, and single photons in pure states a 
resource for others. While great progress has been made 
in the generation of single photons on-demand^], the 
process of parametric downconversion (PDC) constitutes 
a practical room-temperature source of non-classical 
radiation exhibiting unparalleled flexibility for both sets 
of applications. 

Important classes of entangled two-photon states in- 
clude those exhibiting positive frequency correlations [2, 
0| , required for certain metrological applications, and 
those characterized by vanishing correlations in all de- 
grees of freedom, crucial for applications that require the 
concatenation of several sources, such as linear optical 
quantum computation^]. In most PDC sources to date, 
however, the photon pair directions and frequencies are 
strongly correlated, leading to space-time entanglement 
in the joint quantum state. This means that when one 
photon is detected, the other is projected into a mixed 
state, which will not interfere with other photons0,0- 
Experimental proposals for generating factorable states 
have been recently reportedjaSIS; these tend to rely on 
favorable material properties occurring at specific wave- 
lengths. A more general method is to engineer the pho- 
tonic state by coherent manipulation of the probabil- 
ity amplitudes for generating photons of particular fre- 
quency and/or direction. 

The synthesis of two-photon states characterized by 
a joint amplitude of arbitrary shape in two-dimensional 



signal and idler frequency space Lu s -uji requires the use 
of a broadband pump. A narrow bandwidth pump cen- 
tered at (jj p with width e gives rise to photon pairs whose 
frequencies are nearly perfectly anti-correlated such that 
\(jj s +u>i— u) p \ < e/2. Thus, for example, broadband states 
exhibiting vanishing or positive spectral correlations can- 
not be constructed within such a one-dimensional sub- 
space; use of a broadband pump frees us from this restric- 
tion. Work by Klyshko and others^, has explored 
the potential for PDC state engineering of compensating 
dispersive effects in nonlinear crystals by appropriately 
chosen linear media. Here we build upon this pioneering 
work, presenting both theoretically and experimentally 
a technique which relies on a "superlattice" of nonlinear 
and anisotropic linear segments, pumped by an ultra- 
short pulse, which permits the generation of photon pairs 
with an arbitrary degree of entanglement via appropri- 
ate group-delay engineering. When used in waveguided 
geometries, where spectral and spatial degrees of free- 
dom are decoupled|ll| making large collection apertures 
possible, the technique presented here can lead to bright 
photon pair sources with customizable modal character. 

The key to this method is an extension of the no- 
tion of group velocity matching that was introduced 
for the generation of spectrally decorrelated two-photon 
states [3, Il2| , to the idea of manipulating overall group 
delay mismatches between the various fields in structured 
materials. Let us note that a pulsed pump provides a 
timing reference for the PDC photons, not present in the 
CW case. Indeed, the ability to arbitrarily select the ori- 
entation in uj s -LOi space of the joint spectral amplitude 
(thus permitting an arbitrary degree of entanglement) 
requires, in the pulsed case, control over two group de- 
lays (pump-signal and pump-idler) compared to a single 
group delay (signal-idler) in the CW case. Therefore, in 
general, two-dimensional group delay mismatch control is 
required, which can be achieved by interspersing nonlin- 
ear superlattice segments with linear segments exhibiting 
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an appropriate birefringence. 

To be concrete, the two-photon state produced by 
spontaneous PDC for fixed propagation directions^ is: 



\9)=A 
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where A is a normalization constant, aliy:^) (with 
/i = s,i) is the creation operator for the two PDC 
modes and f{uj s , 0Ji) is the joint spectral amplitude (JSA) 
given by the product of the pump envelope function 
(PEF) ajujs+oji) and the phasematching function (PMF) 
cj)(uj s ,uji) []Jj |. The normalized joint spectral intensity 
(JSI) \f(uj s ,LOi)\ 2 represents a probability density func- 
tion for PDC emission at specific frequencies u> s and u>i. 

Our source design is based on a periodic assembly of 
standard crystals and linear birefringent spacers. 
Consider the PDC source shown in Fig. 2] consisting of 
N identical x crystals of length L (cut and oriented 
for degenerate collinear type-II PDC centered at ojq) and 
N — 1 identical linear spacer crystals of length h. 
The PMF for such a superlattice is given by [E0 : 

M»; »i) « Sin[ ^ S, f/i ] sinc [LAk{ Vsi Vi )/2] (2) 

where v s _i — to Si i ~ loq, 

$(i/ s , = LAk(v s ,Vi) + hAn(v s , ^), (3) 

and where Ak{v s ,Vi) and AK,(y Sl Vi) are the phase mis- 
matches in the crystal and spacer segments: 

Ak(v s ,Vi) = k p (v s + Vi) - k s (v s ) - ki(vi) 

AK,(v s ,Vi) = K p (v s + Vi) - K s (v s ) - Ki(vi). (4) 




FIG. 1: Nonlinear crystal superlattice, pumped by an ultra- 
short pulse, yielding an arbitrary degree of entanglement. 

Hence, the total PMF is composed of the product of 
two separate functions: one due to each individual crystal 
segment and one due to the superlattice. A change of 
variables to v± — 2~^(y s ± i>i) and a Taylor expansion 
of the crystal and spacer phasemismatches yields: 

vJ) = /iAk (0) + LAfc (0) + t+v+ + t-V-, (5) 
where Ak^ and Ak' ' are evaluated at ujq and: 

r+ = 2-3 [L(k' s + k' i - 2k' p ) + h{n' s + K'i - 2k' p )] 
r_ = 2-3 [L(k' s - k[) + h{n' s - k[)] . (6) 



Here k'^ and k' (/i = p, s, i) denote the reciprocal group 
velocities (RGV) in the spacer and nonlinear material 
respectively. We assume that Ak^ = 0, thus ensuring 
phasematching in the crystal segments. The resulting su- 
perlattice JSA can be engineered by noting that whereas 
the orientation of the single crystal contribution on the 
v s -Vi plane is fixed, that of the crystal sequence can be 
adjusted by varying the relative thicknesses of the crystal 
and spacer segments 6] . Our approach is to let the single 
crystal PMF be much broader than the spectral structure 
imposed by the superlattice, so that the latter determines 
the resulting two-photon state properties. Under these 
circumstances and modelling the PEF a(uj s + (Ji) as a 
Gaussian function with width a we obtain: 



MNHu^ )/2] 
sin[$(z/+,z/_)/2] PL +/ J V ' 



with $(^+1 v ) given by Eq.El Let us note that rotating 
by 180° subsequent crystal segments, thus flipping the 
incidence angle relative to the optic axis, it is possible to 
limit the maximum transverse walkoff to its single cyrstal 
value (spacer segments exhibit no walkoff for propaga- 
tion parallel or orthogonal to the optic axis). Alterna- 
tively, a waveguided geometry evidently suppresses spa- 
tial walkoff. Simultaneous temporal and spatial matching 
between pump, signal and idler is indeed a powerful tool 
for the engineering of two-photon states. 

The JSA given by Eq. may be adapted to produce 
various photonic states. This may be illustrated by an 
example. States with symmetric JSAs [i.e. f(u s ,u^i) = 
f(u>i,u> s )] exhibit perfect visibility in two-photon Hong- 
Ou-Mandel (HOM) interference, an essential ingredient 
in many optical QIP protocols. Eqns. [7| and [S] reveal 
that by imposing the condition r_ = 0, v_ dependence is 
suppressed; the former leads to a symmetric JSA, in par- 
ticular one exhibiting frequency anticorrelation (which 
may be very tight for a long superlattice) which satisfies 
the HOM symmetry condition despite the fact that, indi- 
vidually, each type-II crystal does not. In our approach, 
the spectral(temporal) distinguishing information from 
individual crystal segments is eliminated by compensat- 
ing the crystal signal-idler temporal walkoff L(k' s — k^) by 
that in the linear birefringent spacers /i(k s — k'A, which 
ensures the suppression of the overall signal-idler tempo- 
ral walkoff term r_. Thus, crystal and spacer materials 
must be chosen with opposite signed signal-idler walkoffs 
while segment thicknesses must be selected appropriately. 

A JSA of this type can be engineered using BBO crys- 
tal segments cut for degenerate collinear type-II PDC at 
800nm (with 42.4° cut angle) and quartz spacers (BBO 
and quartz fast axes oriented orthogonally). For BBO 
segments k' e — k' Q = — 194fs/mm, while for quartz spacers 
K 'e ™ K 'o — 31fs/mm; thus, the spacer to crystal thick- 
ness ratio must be h/L w 6.1 in order that t_ = (the 
plots in Fig. [3 assume L — 250/im and h = 1.53mm). 
Fig. |2I A) shows the single crystal PMF, which does 
not satisfy the HOM symmetry condition jl4|. Fig.[2lB) 
shows the superlattice PMF contribution (left-hand fac- 
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FIG. 2: Density plots of: (A) single crystal segment PMF exhibiting lack of symmetry, (B) BBO-quartz-BBO superlattice 
contribution to the PMF, (C) pump envelope function, (D) JSI symmetrized by the superlattice. (E)Time (with respect to 
pump pulse) at which signal/idler generated halfway through the first crystal traverse each point of two crystal segments. 
(F)Same as before for a BBO-quartz-BBO superlattice; signal-idler walkoff is limited to its single crystal value. 



tor in Eq. |2J). Fig. |2{C) shows the PEF (centered at 
400nm with 1.95nm bandwidth). Fig.^D) shows the re- 
sulting crystal-spacer-crystal superlattice JSI exhibiting 
the desired symmetry; the single crystal spectral band- 
width is sufficiently large that the superlattice determines 
the photon pair properties. 

The group delay engineering above may be understood 
by considering how the quartz spacer limits the signal- 
idler temporal walkoff. Fig. |2[E) shows for each point 
along a single 500/^m thick crystal (considered as two 
250/^m crystals in sequence) the time, with respect to 
the pump pulse, at which each of the signal and idler 
wavepackets traverses each point of the crystal, assuming 
a photon pair created halfway through the first crystal. 
Fig. EfF) shows the effect of the quartz spacer: it limits 
the maximum signal-idler walkoff to that observed in a 
single 250/^m crystal (it increases, however, the pump- 
PDC walkoff, which does not impact the HOM symme- 
try condition). Thus, even though signal-idler temporal 
walkoff cannot be entirely suppressed, it can be limited 
in a superlattice to an arbitrarily small value. 

In order to illustrate the power of this concept, we im- 
plemented the simple superlattice described above. The 
experimental configuration used [see Fig. [3] relics on a 
collinear type II HOM interferometer (HOMI) . A train of 
ultrashort pulses centered at 405nm (2nm FWHM band- 
width, 82MHz repetition rate and 150mW power), ob- 
tained by frequency doubling a mode-locked Ti:Sapphirc 
beam, was weakly focused onto a collinear type II PDC 
source. The latter is based on 250imr length BBO crys- 
tal segments, and could be configured as single-crystal, 
double-crystal or as a crystal-quartz-crystal superlattice. 
A polarizing beamsplitter (PBS1) established a Michel- 
son interferometer, each polarization in one arm. Quar- 
ter wave plates (QWP1/QWP2) reverse the polariza- 
tion, so that both photons exit PBS1 from its second 
port. One of the arm lengths was adjustable via a mo- 
torized translation stage (TS). The PDC exiting the in- 
terferometer was coupled into a single mode fiber (SMF) 
which ensured a high degree of signal-idler spatial mode- 
matching, prior to which a dichroic beamsplitter (blue re- 
flecting/red transmitting) followed by a red-colored glass 
filter, separated the PDC light from the blue pump. 
A half-wave plate (HWP) oriented the orthogonal sig- 
nal/idler modes at 45° with respect to the axes of a sec- 
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FIG. 3: BB01/BB02: 250rtm BBO crystals, QTZ: 1.53mm 
quartz spacer, PBS1/PBS2: polarizing beam splitters, 
QWP1/QWP2: quarter waveplates, TS: translatable mirror, 
MO: 16x microscope objetive, SMF: single mode fiber, HWP: 
half wave plate, BP1/BP2: bandpass filters of 40nm width, 
APD1/APD2: avalanche photo diodes. 



ond polarizing beam splitter (PBS2) which formed the 
basis for a collinear HOMI. Each output mode from the 
PBS was focused on the active area of a silicon photon- 
counting module; wide-bandpass filters (wider than the 
PDC bandwidth) were used for background suppression. 
Data was taken by monitoring coincidence counts as a 
function of the translation stage (TS) position. 

The data in Fig. 0] shows coincidence count rate data, 
normalized to a unit background, as a function of the 
signal-idler relative delay. Solid curves represent theo- 
retical curves obtained from a calculation similar to that 
in Ref. Fig. UfA) shows the interference curve re- 

sulting from a single crystal segment. Fig. 0IB) shows 
a broadened dip with reduced visibility, due to the in- 
creased signal-idler temporal walkoff [quantified by r_, 
Eq. resulting from two crystal segments in sequence. 
Fig. E^C) shows the result of placing a 1.6mm length 
quartz spacer between the crystals [see Fig. [3]. Note that 
the visibility and dip width are restored to their single 
crystal values, thus demonstrating the effect of limiting 
the maximum walkoff to that experienced by a single 
crystal. Fig.0jD) shows the effect of rotating the quartz 
spacer by 90°, thus reversing the extraordinary and or- 
dinary character in the spacer. In this latter case, the 
superlattice contribution is rotated away from the diag- 
onal axis on the u) s — a;, plane, thus destroying fulfilment 
of the HOM symmetry condition; the interference curve 



4 



1.0 











OJ 


" ' \ * / * " 


\ /" 
\ 
*\ • 


coincidence ra- 


\ / 

(B) 


-. •, . 








1.0 

01 

2o.8 

OJ 

s a6 

xs 




•\ /• 


'V # •V* 




£0.4 




V* (0 


'o 

^0.2 


(D) 



-0.05 0D5 0.1 -0.05 0.05 0.1 



delay [ps] delay [ps] 

FIG. 4: HOM interference curves, normalized to a unit back- 
ground to suppress source brightness variation resulting from 
re-alignment for each experiment. (A) single crystal segment. 
(B) Two crystal segments in sequence, exhibiting reduced visi- 
bility and broadening. (C) BBO-quartz-BBO superlattice ex- 
hibiting restoration of visibility. (D) BBO-quartz-BBO with 
quartz rotated by 90°; visibility no longer restored. 

now contains two dips, each associated with one of the 
crystals. Let us note that our conclusions depend on an 
analysis of HOM interference features and are indepen- 
dent of source brightness; however, a future refinement 
is the careful selection of the spacer thickness in order to 
rule out thickness-dependent count rate oscillations 

So far we have considered the effect of signal-idler tem- 
poral walkoff. In general the pump experiences a different 
group velocity from the signal and idler photons, so that 
there is a corresponding temporal walkoff between the 
pump and the PDC photon pairs. It is evident that the 
condition r + = (see Eq. |SJ) closely resembles the group 
velocity matching (GVM) condition|3, ^3 by which the 
pump RGV is matched to the mean RGV of the sig- 
nal and idler wavepackets. In a superlattice, however, 
it is the total group delay mismatch which must van- 
ish. This generalized condition can lead to spectrally 
factorable states, or to states exhibiting positive spectral 
correlations. For example, r + = can be satisfied us- 



ing type-II BBO crystal segments and calcite spacers |6J. 
Note that the GVM condition fulfilled by the superlattice 
depends crucially on L and h. The thickness error toler- 
ance becomes more restrictive for longer assemblies (i.e. 
for higher N); for example for a superlattice of 10 BBO 
crystals and 9 calcite spacers designed to yield a positive 
correlation state, the tolerance is about 2%. In addition, 
the technique presented here can be adapted to produce 
states exhibiting a very large Schmidt number (and thus 
a high degree of space-time entanglement) |l5j which con- 
stitute extreme examples of non-classical correlations [l6j . 

In conclusion, we have shown that a periodic ar- 
ray of nonlinear crystals and linear birefringent spacers, 
pumped by an ultrashort pump, can result in PDC pho- 
ton pairs exhibiting an engineered JSA function. We 
have illustrated this concept by demonstrating a source 
capable of high-visibility HOM interference, in a type- 
II pulsed pump configuration, using using materials that 
do not generate such states naturally. Our scheme has 
wider applicability; it can generate a broad class of two 
photon states including those with factorable and posi- 
tively correlated JSAs, as well as those with a very large 
Schmidt number. It is conceivable that appropriate su- 
perlattices could be implemented in monolithic nonlinear 
optical structures where locally- modified dispersion char- 
acteristics are obtained by ion-exchange, doping or other 
mechanisms. This would open up a powerful new source 
engineering technique, akin to quasi phasematching but 
in the group velocity, rather than phase velocity, domain. 
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